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Stars with short orbital periods at the center of our galaxy offer a powerful 
and unique probe of a supermassive black hole. Over the past 17 years, the 
W. M. Keck Observatory has been used to image the Galactic center at the 
highest angular resolution possible today. By adding to this data set and ad- 
vancing methodologies, we have detected SO-102, a star orbiting our galaxy's 
supermassive black hole with a period of just 11.5 years. SO-102 doubles the 
number of stars with full phase coverage and periods less than 20 years. It 
thereby provides the opportunity with future measurements to resolve degen- 
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eracies in the parameters describing the central gravitational potential and to 
test Einstein's theory of General Relativity in an unexplored regime. 

The recent advent of high resolution imaging capabilities on large ground-based telescopes 
has afforded us a view of the center of our galaxy that is substantially altering our understanding 
of this region and by analogy the centers of other galaxies. The earliest high-resolution images 
in the near-infrared (NIR) revealed the presence of a black hole with a mass 4 • 10 6 times that 
of the sun by detecting stars at the very heart of the galaxy, the so called S-stars, with velocities 
of up to ^10,000 km/s (7, 2). The detection of accelerations in the motion of these stars (3, 4) 
further strengthened the case for a non-stellar dark mass. After passing periapse - the closest 
point to the black hole in the orbit of a star - a Keplerian orbit could be determined for the star 
SO-2 (with an orbital period of 16 years), proving the black hole hypothesis beyond reasonable 
doubt (5, 6). This star can now be traced in its motion around the black hole with such precision 
and complete phase coverage that its three-dimensional orbit gives the most accurate description 
to date of the central mass and distance from our Solar system to the center of the Galaxy (7-9); 
no other star has so far been reported with more than ^40% of its orbit covered by data. 

Here, we report the detection of another star, SO- 102, in orbit around the supermassive black 
hole at the center of our galaxy (70). It has an orbital period of 1 1.5 years, which makes it the 
shortest period star known yet. 

The data sets that we used for this study come from high-resolution imaging observations of 
the Galactic center that were obtained at the W. M. Keck Observatory between 1995 and 2012 
with speckle imaging and adaptive optics. Speckle imaging with the near-infrared imaging 
camera NIRC (77, 72) on Keck I provided the earliest high angular resolution imaging data sets 
(1995 - 2005), which are reported with their observational setup in several earlier papers (2, 3, 
7, 13, 14). In summary, several thousand data frames capture images of the Galactic center on 
timescales of ~0.1 sec, which is short enough compared to the coherence time of atmospheric 
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turbulence to freeze the distorting effects of the Earth's atmosphere and preserve the source's 
information at all spatial frequencies transmitted by the telescope. Table SI summarizes key 
observational information for the 26 epochs of speckle imaging data that are available. 

Adaptive optics (AO) observations of the Galactic center were carried out between 2004 and 
2012, with substantially more data per epoch being obtained in 2006 and onwards. Previous 
publications contain the observational details for the first 9 epochs of AO observations (7, 14- 
16). We collected 12 further epochs of AO observations between 2008 May and 2012 May, 
using an identical observational set-up to all AO observations obtained in 2006 and onwards. 
For all but two epochs, the high-order AO compensation was based on observations of a bright 
(R^IO mag), on-axis laser guide star (LGS), generated with the Keck II sodium laser, and 
the low-order (e.g. tip-tilt) AO corrections were determined from measurements of the natural 
guide star (NGS) USNO 0600-28577051 (R = 13.7 mag and Ar Sgr A* = 19"). Table S2 provides 
the details of the 20 AO epochs that are used here. 

The data analysis for this study falls into several distinct stages (77). First, for each epoch 
of observations, we constructed an average image of all the high quality data collected as well 
as three sub-set images, each composed of one third of the data (see Fig. 1 for a zoom into 
the region of an AO image where SO- 102, SO-2, and the electromagnetic manifestation of the 
black hole, Sgr A*, reside). Second, both previously known and newly detected stars were 
identified in the images, along with their photometric and astrometric characteristics, with the 
point spread function (PSF) fitting program StarFinder (18). Third, the positions of the stars in 
each epoch were transformed into a common coordinate system. Fourth, orbital model fits to 
the positions of stars with significant orbital phase coverage provided estimates of each star's 
orbital parameters. 

The gravitational potential in this region of the Galaxy is dominated by the black hole, 
and to a good approximation the stars follow Keplerian orbits. While deviations from a pure 
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Keplerian orbit will likely be detected in the future, this cannot be expected with the current 
data set, and we therefore assumed a point-mass potential and Newtonian gravity. Within this 
framework, we used a x 2 -minimization routine to find the best fit orbital parameters (7). The 
six orbital elements describing the star's trajectory are the period, eccentricity, time of closest 
approach (periapse passage), inclination, angle to periapse, and the angle of the ascending node. 
The seven parameters describing the gravitational potential are mass, 3-dimensional position 
(distance from Earth & focal point of the orbital ellipse) and 3-dimensional velocity of the 
point-mass. 

We have covered a complete orbit of SO- 102 with astrometric observations. Therefore, the 
data contain enough information that they can in principle be used to fit for both the star's orbital 
elements and also for the parameters of the potential (79). However, this requires the same low 
level of systematic effects in the data as for SO-2, which dominates our current knowledge 
about the potential and is a factor of 16 brighter than SO- 102. Considering that fainter sources 
are much more prone to source confusion than brighter ones, it is unlikely that SO- 102 has the 
same low level of systematic errors (7, 8). 

A way of minimizing the effects of possible confusion is to not fit for the potential, but rather 
fix it to the values that have been derived from the orbit of SO-2 (Table S4). This left only S0- 
102's orbital elements free, which were then found by fitting the data (Fig. 2). Overall, the data 
are well fit by this model, which has a reduced-^ 2 of 2.0. The most discrepant points are the two 
epochs in 2009, both of which lie 2.3a away from the best fit. This is the region around closest 
approach to the black hole, where Sgr A* and SO- 104 (20), which are blended in this epoch, 
are very close (^55 mas) to SO- 102. This constellation most likely leads to an astrometric bias, 
which supports our conservative approach of fixing the parameters of the potential. 

We determined uncertainties on S0-102's orbital parameters via Monte Carlo simulations 
assuming Gaussian error statistics. Using the best fit and the error bars on the measured data 
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points, we generated 10 5 artificial data sets and fitted an orbit to each realization. The proba- 
bility distribution functions for the orbital elements are the distribution of best fit values from 
the Monte Carlo simulations (Table 1). This approach assumes that the positional errors are sta- 
tistical in nature and does not account for possible systematic contributions like unrecognized 
source confusion. 

Our high angular resolution imaging campaign of the Galactic center started in 1995, con- 
stituting a time baseline of 17 years, 7 of which were carried out with deep AO observations. 
While this is long enough to detect accelerations (i.e. curvature) of several stars within the 
central arcsecond, it only allows a reliable orbit determination for the shortest period stars. We 
consider an orbit reliable if it is likely that an increase in data and a decrease in noise does not 
lead to major changes in the orbital elements. As a rule of thumb, 50% of the orbit needs to be 
covered by observations in order to determine the orbit reliably (27). Our data set traces SO- 102 
through a complete orbit, making it the only star other than SO-2 for which such a fraction of 
the orbit is sampled. The key property that enabled this orbital coverage is S0-102's period. 

With an orbital period of 11.5 years, SO- 102 has the shortest period among all known stars 
orbiting the black hole in our galaxy's center (Fig. 3). It is 5 years (~30%) shorter than S0-2's 
orbital period, which was the previously known shortest-period star. The period (or, alterna- 
tively, the semi-major axis, which is related to the period via Kepler's law) is the most important 
orbital parameter because it not only makes it possible to sample a substantial fraction of the 
orbit by observations, but also simplifies the detection of relativistic effects, most of which are 
cumulative and build up with increasing phase coverage until they breach the detection thresh- 
old. 

A test of Einstein's General Theory of Relativity is the next goal in Galactic center research 
now that the existence of a black hole is well established. This theory has so far passed all tests 
on solar system scales with flying colors. However, the gravitational field of the sun is very 
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weak, with typical experiments probing regimes of up to e ~ GMj (Re 2 ) ~ 10~ 6 . Here, G is 
the gravitational constant, M the mass scale, R the distance scale, and c the speed of light (22). 
The gravitational fields that have been probed in tests using double neutron stars like the Hulse- 
Taylor binary pulsar are of the same magnitude, since the masses and separation of the neutron 
stars are comparable to the mass and radius of the Sun (23). In the Galactic center, stars like 
SO- 102 and SO-2 probe gravity regimes that are two orders of magnitude stronger, e ~ 10~ 4 . 

The effects of curved space time manifest themselves in two ways: the orbit of a star deviates 
from its Keplerian approximation, and the path and wavelength of the light emitted by it get 
changed. It has been shown that for SO-2 the relativistic redshift should become measurable 
at S0-2's next periapse passage in 2018 (24, 25). Although this is a measurement of S0-2's 
redshift, SO- 102 will render this measurement more precise. The reason for this is that, with AO 
observations of S0-102's apoapse passage in the next years, there will be enough information in 
the data set to lower the level of systematic effects to an insignificant level, and thereby allow 
SO- 102 to be used to constrain the central gravitational potential. This additional information 
will then translate into a more precise determination of S0-2's orbital elements and therefore 
redshift measurement. 

The relativistic redshift will also become observable in SO- 102 itself. For SO- 102, the dif- 
ference between the purely Newtonian and fully relativistic Doppler shifts amounts up to ^90 
km/s. This is within the capabilities of upgraded or next-generation near-infrared spectrographs 
for sources of this magnitude. Since the Keplerian orbital parameters must be known precisely 
to infer the relativistic contribution to the redshift, only stars with well sampled astrometric 
orbits can be used to detect deviations from a Newtonian orbit. Additionally, the relativistic 
redshift is strongest at periapse passage, and SO-2 and SO- 102 are the only known stars that 
pass periapse within the next 10 years and that will have an orbital phase coverage greater than 
50% then. Since SO- 102 will pass closest approach 3 years after SO-2, an observation of its 
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relativistic redshift will provide a crucial check on systematic effects in S0-2's result, and test 
the independence of the gravitational redshift from the interior structure of the star - a key 
assumption of General Relativity. 

Other than the gravitational redshift, the leading order General Relativistic effect is the 
precession of the periapse within the orbital plane (26-31). This is a cumulative effect that adds 
up from orbit to orbit and is therefore much more likely to be detected in stars that complete 
several orbits during the timescale of an observational campaign. Angelil et al. (30) show that 
the effect of periapse precession scales as period' 1 , quantifying the significance of a short 
period for the test of General Relativity in an unprecedented regime. 

The periapse precession, which leads to rosetta-shaped orbits, consists of two components: 
a General Relativistic part that leads to a prograde precession, and a Newtonian part that leads 
to a retrograde precession. The Newtonian deviation of a Keplerian orbit is caused by the 
inevitable presence of extended dark mass around the black hole, which means that the grav- 
itational potential is not point-like (27, 31, 32). This extended mass distribution is thought to 
mainly consist of stellar remnants - stellar mass black holes and neutron stars - supplied by the 
dense nuclear star cluster (33, 34). In order to measure the General Relativistic part only, the 
detection of the periapse shift in more than one star is required to break the degeneracy. This 
makes a star like SO- 102 a necessary ingredient for a measurement of the warping of spacetime 
around a supermassive black hole. 

Given the astrometric precision with which a star can be located on the detector (^0.1 mas 
and ~1 mas for SO-2 and SO- 102, respectively) and with which a stable absolute reference 
frame can be constructed (currently ~0.1 mas/yr), the precession of both S0-2's and S0-102's 
apoapse - which amounts to ^1 mas and ^0.3 mas, respectively - is likely only detectable with 
the next generation of telescopes like the Thirty Meter Telescope. In this upcoming era, SO- 102 
will play as dominating a role as SO-2 plays in today's Keck era. 



7 



References and Notes 

1. A. Eckart, R. Genzel, Mon. Not. R. Astron. Soc. 284, 576 (1997). 

2. A. M. Ghez, B. L. Klein, M. Morris, E. E. Becklin, Astrophys. J. 509, 678 (1998). 

3. A. M. Ghez, M. Morris, E. E. Becklin, A. Tanner, T. Kremenek, Nature 407, 349 (2000). 

4. A. Eckart, R. Genzel, T. Ott, T., R. Schodel, Mon. Not. R. Astron. Soc. 331, 917 (2002). 

5. R. Schodel et al, Nature 419, 694 (2002). 

6. A. M. Ghez et al, Astrophys. J. 586, L127 (2003). 

7. A. M. Ghez et al, Astrophys. J. 689, 1044 (2008). 

8. S. Gillessen et al, Astrophys. J. 692, 1075 (2009). 

9. S. Gillessen et al, Astrophys. J. 707, LI 14 (2009). 

10. SO- 102 might be identical to the star labeled S55 in (8), for which only linear motion is 
reported there. 

11. K. Matthews, B. T. Soifer, in Astronomy with Infrared Arrays: The Next Generation, ed. 
I. McLean (Kluwer, 1994), Astrophysics and Space Science, v. 190, p. 239 

12. K. Matthews, A. M. Ghez, A. J. Weinberger, G. Neugebauer, The Publications of the 
Astronomical Society of the Pacific 108, 615 (1996). 

13. A. M. Ghez et al, Astrophys. J. 620, 744 (2005). 

14. J. R. Lu et al, Astrophys. J. 690, 1463 (2009). 

15. A. M. Ghez et al, Astrophys. J. 635, 1087 (2005). 

8 



16. S. D. Hornstein et al, Astrophys. J. 667, 900 (2007). 

17. A detailed description of the methods is available as supplementary material on Science 
Online. 

18. E. Diolaiti, et al., Astron. Astrophys. Suppl Sen 147, 335 (2000). 

19. The distance to the black hole and its line-of-sight velocity cannot be determined, because 
a measurement of the star's radial velocity is needed for this. SO- 102 is currently just 
beyond the spectroscopic limit. 

20. SO- 104 is a previously unknown detected star; its orbital motion is presented in the sup- 
plementary material on Science Online. 

21. W. I. Hartkopf, B. D. Mason, E. Worley, Astron. J. 122, 3471 (2001). 

22. C. Will, Living Reviews in Relativity (available at http://www.livingreviews.org/lrr-2006- 
3), 2001 

23. The strength of the gravitational field we quote here relates to the inter-body potential 
and therefore orbital velocities. However, we want to note that the strong-field internal 
gravity of the neutron stars could leave imprints on the orbital motion. 

24. S. Zucker, T. Alexander, S. Gillessen, F. Eisenhauer, R. Genzel, Astrophys. J. 639, L21 
(2006). 

25. R. Anglil, P. Saha, Astrophys. J. 734, L19 (2011). 

26. M. Jaroszynski, Acta Astronomica 48, 653 (1998). 

27. G. F. Rubilar, A. Eckart, Astron. Astrophys. 374, 95 (2001). 



9 



28. N. N. Weinberg, M. Milosavljevic, A. M. Ghez, Astrophys. J. 622, 878 (2005). 

29. C. M. Will, Astrophys. J. 674, L25 (2008). 

30. R. Anglil, P. Saha, D. Merritt, Astrophys. J. 720, 1303 (2010). 

31. D. Merritt, T. Alexander, S. Mikkola, C. M. Will, Phys. Rev. D 81, 062002 (2010). 

32. A. Gualandris, S. Gillessen, D. Merritt, Mon. Not. R. Astron. Soc. 409, 1146 (2010). 

33. M. R. Morris, Astrophys. J. 408, 496 (1993). 

34. J. Miralda-Escude, A. Gould, Astrophys. J. 545, 847 (2000). 

35. S. Yelda et al, Astrophys. J. 725, 331 (2010). 

36. W. I. Clarkson et al, Astrophys. J. 751, 132 (2012). 

37. S. D. Hornstein et al, Astrophys. J. 577, L9 (2002). 

38. S. D. Hornstein, thesis, UCLA (2007). 



39. R. Schodel et al, (available at http://arxiv.org/abs/! 1 10.2261v3 ) 



40. J. C. Christou, S. T. Ridgway, D. F. Buscher, C. A. Haniff, D. W. McCarthy Jr., in Astro- 
physics with Infrared Arrays, ASP Conference Series, Volume 14, 133 (1991). 

41. J. Primot, G. Rousset, J. C. Fontanella, Journal of the Optical Society of America A 1, 
1598 (1990). 

42. M. G. Petr, V. Coude Du Foresto, S. V. W. Beckwith, A. Richichi, M. J. McCaughrean, 
Astrophys. J. 500, 825 (1998). 

Acknowledgements: We thank the staff of the Keck Observatory for all their help in ob- 
taining the new observations. Support for this work was provided by NSF grant AST 

10 



0909218, the Levine-Leichtman Family Foundation, and the W. M. Keck Foundation. 
RS acknowledges support by the Ramon y Cajal programme, by grants AYA201 0-1 7631, 
AYA2009- 13036, and PA1002584 of the Spanish Ministerio de Economia y Competitivi- 
dad, and by grant P08-TIC-4075 of the Junta de Andalucia. The W. M. Keck Observatory 
is operated as a scientific partnership among the California Institute of Technology, the 
University of California and the National Aeronautics and Space Administration. The 
Observatory was made possible by the generous financial support of the W. M. Keck 
Foundation. The authors wish to recognize and acknowledge the very significant cultural 
role and reverence that the summit of Mauna Kea has always had within the indigenous 
Hawaiian community. We are most fortunate to have the opportunity to conduct observa- 
tions from this mountain. The data described in the paper are presented in the Supporting 
Online Material. 



11 



Table 1 . Orbital elements for SO- 1 02 a 



Parameter [Unit] 


Value 


SO-102's orbital parameters 


Period [years] 


11.5 ±0.3 


Time of Closest Approach [year] 


2009.5 ±0.3 


Eccentricity 


0.68 ±0.02 


Inclination 13 > c [degrees] 


151 ±3 


Angle to Periapse [degrees] 


185 ±9 


Position Angle of the Ascending Node c [degrees] 


175 ±5 


Parameters of the potential 


Mass [10 6 M Sun ] 


4.1 ±0.4 


Distance [kpc] 


7.7 ±0.4 



a The best fit has a x 2 of 39.96 with 20 degrees-of-freedom. 
b 90 deg is edge-on, and deg is face-on. 

c The allowed ranges for inclination and angle of the ascending node 
are [0,180]. 

d The parameters that describe the gravitational potential have been 
taken from S0-2's orbit. We list here mass and distance, see Table S4 
for 2D position and 3D velocity of the central mass. 
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Figure 1 A Keck/NIRC2 adaptive optics image from May 2010 showing the short-period star 
SO- 102, which is besides SO-2 the only star with full orbital phase coverage, and the electro- 
magnetic counterpart of the black hole, Sgr A*. The image was taken at a wavelength of 2.12 
/xm and shows the challenge of detecting SO- 102, which is 16 times fainter than SO-2 and lies 
in this crowded region. 
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Orbits of SO-2 & SO-102 
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Figure 2 The orbits of SO-2 (black) and SO-102 (red). The data points and the best fits are 
shown. Both stars orbit clockwise. The dashed lines represent the part of the orbits that have 
been observed with Speckle data, the solid lines indicate adaptive optics observations. The 
data points for SO-2 range from the year 1995 to 2012, S0-102's detections range from 2000 to 
2012. The connecting lines to the best fit visualize the residuals. Note that while the best-fit 
orbits are not closing, the statistically allowed sets of orbital trajectories are consistent with a 
closed orbit. SO-102 has an orbital period of 11.5 years, 30% shorter than SO-2, the previously 
known shortest-period star. 
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Figure 3 The orbital periods and magnitudes of all known stars orbiting the Galactic black hole 
(7, 8, this publication). The size of the points is scaled with the orbital phase that is covered 
by observations. The higher the fraction of the orbit that is sampled by observations the more 
reliable the orbital solution is. The three vertical lines indicate, from left to right, the limiting 
magnitude for Speckle shift-and-add (SAA), Speckle holography (SpH), and adaptive optics 
(AO) data sets. We define the limiting magnitude as the median of the limiting magnitudes 
of individual epochs, which are listed in Tables S2 & S3. The filled point represents the star 
SO-102. 
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